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Abstract 

We report on a lattice QCD calculation of the irNa term, the scalar form 
factor, and (iV|ss|iV). The disconnected insertion part of a n N is found to 
be 1.8 db 0.1 times larger than the connected insertion contribution. The q 2 
dependence of a 7T N(q 2 ) is about the same as GE{q 2 ) of the proton so that 
(T nN {2ml )-a nN {0) = 6.6±0.6MeV. The ratio y = (N\ss\N) / (N\uu+ dd\N) = 
0.36 ± 0.03. Both results favor a ct^n ~ 53 MeV, slightly larger than our direct 
calculation of a n N = 49.7 ± 2.6 MeV. We also compute F s and D s and find 
that the agreement with those from the octect baryon mass splittings crucially 
depends on the inclusion of the large disconnected insertion. Finally, we give 
our result for the KNa term. 
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Like the pion mass in the meson sector, the irNa term is a measure of the explicit 
chiral symmetry breaking in the baryon sector. It is considered a fundamental quan- 
tity which pertains to a wide range of issues in the low-energy hadron physics, such 
as quark and baryon masses, strangeness content of the nucleon, pattern of SU(3) 
breaking, irN and KN scatterings, kaon condensate in dense matter, trace anomaly, 
and decoupling of heavy quarks. Defined as the double-commutator of the isovec- 
tor axial charge with the hamiltonian density taken between the nucleon states, i.e. 
a wN = l/3£o=i,3(AT|[Qo> [Q 5 a ,H(0)}}\N) which appears in the off-shell ttN scattering 
amplitude [[lj], it has in QCD the expression 

a wN = m(N\uu + dd\N), (1) 

where m = (m u + m<j)/2. 

It is shown that at lowest order in m n (i.e. m 2 ), it is equal to the unphysical, but 
on-shell, isospin even ttN scattering amplitude at the Cheng-Dashen point, E^ = 
f 2 T + (s = m 2 N ,t — q 2 — 2m 2 ). Thus E^ can be extracted from ttN scattering 
experiment via fixed-t dispersion relation for instance Q . It is further shown || that 
the next higher order term which is nonanalytic in quark mass (i.e. proportional to 
m 3 / 2 or m 3 ) drops out if E^at is identified with a 7r Ar(2m 2 r ) || which is only a function 
of q 2 . This shows that the difference A# in the relation T, n ^ = cr 7r Ar(2m 2 r ) + A^ is 
of the order m^/m% and has been shown to be indeed negligible (~ 0.35 MeV) in a 
chiral perturbation calculation |3|. || . 

Various estimates of E^at have ranged from 22 to 110 MeV over the years, but 
eventually settled around 60 MeV Q. On the other hand, a puzzle was raised by 
Cheng |J. If one assumes that (iV|ss|iV) = 0, a reasonable assumption from the OZI 
rule, the cr^jy obtained from the octet baryon masses gives only 32 MeV, almost a fac- 
tor two smaller than E^ extracted from the 7rN scattering. This puzzle was tackled 
from both ends. First, the scalar form factor was calculated [SI in chiral perturbation 
theory (x PT) with the two correlated pions as the dominating intermediate state. As 
a result, the scalar form factor is found to be exceedingly soft which leads to a large 
change of cr^N^q 2 ) in a small range of q 2 , i.e. Act^at = o'nN(2m 2 [ ) — <J 7T ^(0) = 15.2 ±0.4 
MeV. Thus, this reduces ovjv 

&nN = E^at — A(T nN (2) 

to ~ 45 MeV. The remaining discrepancy between a n N and is reconciled if one 
admits the possibility of a large ss content in the nucleon || ||. From the pattern of 
SU(3) breaking in the octet baryon masses, one finds 0. [5| 

<r vN = a^f/{l-y), (3) 

where y = 2(N\ss\N) / (N\uu + dd\N). Given af ] N = 32MeV from the octet baryon 
masses ||, or 35(5) MeV from the one loop x PT calculation || and a^N = 45MeV, 
eq. (Q) implies y = 0.2 — 0.3. 

Hence, a consistent solution seems to have emerged which suggests that a^N ~ 
45 MeV, Aa n N ~ 15 MeV, and y ~ 0.2 — 0.3. In this letter, we undertake a lattice 



1 



QCD calculation of the above quantities to scrutinize the viability of this resolution. 
It turns out that our study strongly points to a significantly different solution as we 
shall show. 

The calculation of a^jq in lattice QCD has been attempted by several groups [0, |§ 
who employed the Feynman-Hellman theorem 



rrS^- = rh(N\uu + dd\N) c i + m(N\uu + dd\N) D / (4) 
dm 

and obtained a n N through the derivative of the nucleon mass. We note that in eq. 
(f|) the connected insertion (C.I.) part comes from the differentiation with respect to 
the valence quark propagator; whereas the disconnected insertion (D.I.) comes from 
the derivative of the fermion determinant. Their contributions to the scalar current 
ipip in the nucleon are shown schematically in Fig. 1. In the quenched approximation 
approach, it is found that a^N obtained from the derivative of the nucleon mass is 
only about 15 — 25 MeV |7j. This is much smaller than the phenomenological value 
of ~ 45 MeV 0. The smallness of i n this case is traced to the fact that the 
nucleon mass in the quenched approximation is calculated with the determinant set 
to a constant, so that its derivative corresponds to the C.I. only (which is verifiable 
by comparing to the direct evaluation of the C.I. || |K]]) and it does not involve 
D.I. which can be substantial. Indeed, when the derivative of Mn is calculated with 
dynamical fermions included, it is found || that the l.h.s. of eq. (^) which now 
includes the D.I. becomes ~ 2 to 3 times larger than the C.I. contribution. This 
implies a large contribution of the D.I.. Since the error on dM^/dm is quite large 
flO| , ||, we decided to calculate the D.I. directly with the help of the Z 2 noise [|TT |. 



Following our calculation of the flavor- singlet g A we calculate the C.I. and D.I. 
of a^N directly in the quenched approximation. In terms of the Feynman-Hellman 
theorem, it would correspond to calculating dM^ / dm by taking the derivative of the 
determinant first before setting it to a constant. 

Lattice calculations of three-point functions have been used to study the EM []j"3"] , 
axial (isovector) fll~4| , pseudoscalar(7r iViV) Jl5j form factors, and the flavor-singlet g° A 
|T2| |. For the scalar current, we use S(x) = 2k/8k c [uu(x) + dd(x)}, where we have 
implemented the mean-field improvement factor 8k c to define the lattice operator 

The C.I. is calculated in the same way as the isovector axial coupling g\ [|i~4 



Numerical details are given in Ref. [Q. Like in Ref. [O, the lattice renormalized 
9s,con = (N\uu + dd\N)^ on has been calculated for k = 0.154,0.152, and 0.148, cor- 
responding to quark masses of about 120, 200, and 370 MeV respectively (the scale 
a^ 1 = 1.74(10)GeV is set by the nucleon mass), and is plotted in Fig. 2(a). The 
calculations were done on a quenched 16 3 x 24 lattice at (3 = 6.0 with 24 gauge con- 
figurations as in the previous cases [14], |P2"| . Due to the fact that the quenched x PT 
calculation exhibits a leading non-analytic behavior of m 3 ^ 2 for the nucleon mass [[17 



we extrapolate gs con to the chiral limit (k c = 0.1568) with the form C + Dm 1 / 2 . This 
is so because g$ con = dM N /dm in the quenched approximation as we alluded to ear- 
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lier in eq. (|]). As a result, we obtain gg ccm = 3.04(9) as shown in Fig. 2(a). The gs 
in the continuum with the MS scheme is related to its lattice counterpart by the rela- 
tion gs = Zsg$, where Z$ is the finite lattice renormalization constant. The one-loop 
calculation gives Zs = 0.995 for (3 = 6.0 [|16|], from which we find gs,con = 3.02 ± 0.09. 
We also computed isovector g s = (N\uu — dd\N) which does not involve the D.I. and 
find it to be 0.63(7). 

Since a n N is renormalization group invariant, the C.I. contribution is <J n N,con = 
™9s,con where m is the lattice quark mass. From and M N , we find m = 5.84(13) 
MeV. Thus, a^con = 17.8(9) MeV which agrees well with previous calculations 
|7|, ||, |(J. The C.I. part of the form factor is obtained by extrapolating gs con (q 2 ) a ^ 
different k to the chiral limit. It is plotted in Fig. 2(b) together with g\{q 2 ), the 
isovector axial form factor. We see that they are almost identical within errors. In 
so far as the concept of meson dominance goes, this reflects in part the fact that the 
isovector scalar meson and Ai are essentially degenerate in the lattice calculation. 
Like in the case of the axial coupling constants we also find that the ratio 



■R-s — 9% I ' gs,con dips below the SU(6) result of 1/3 as the quark mass becomes lighter. 
This is interpreted as due to the cloud quark/anti-quark effect and is responsible for 



the u - d parton difference reflected in the Gottfried sum rule plf . Only when the 
cloud degree of freedom is eliminated in the valence approximation J2lj where the 
Fock space is limited to the valence do we recover the SU(6) limit. This indirectly 
shows the effect of the cloud quarks in the C.I.. 

We calculate the D.I. in Fig. 1(b) the same way we did for the D.I. part of g\ 
|12|1 by summing over t. For tf >> a, this sum becomes 

r G ™p(tf,P,t,q)Gr P {t,P) _ _ ef ^^ 



t 



G a P a P (tfiP)G a P a P (t,q) 



-* const + t f g Sdia (q ) (5) 



tf»a 



Thus, we calculate the sum as a function of tf and take the slope to obtain the 
D.I. part of gg. Since the D.I. involves quark loops which entail the calculation of 
traces of the inverse quark matrices, we use the proven efficient algorithm to estimate 



these traces stochastically with the Z 2 noise |Tlj which was applied to the study of 
9°A 0- 

The results of eq. © with 300 complex Z 2 noise and 50 gauge configurations 
for k = 0.148,0.152 and 0.154 are presented in Fig. 3. The corresponding gs dis = 
(N\uu + dd\N) dis are obtained from fitting the slopes in the region tf > 8 where the 
nucleon is isolated from its excited states with the correlation among the time slices 
taken into account [0. The resultant fits covering the ranges oft/ with the minimum 
X 2 are plotted in Fig. 3. Finally, the errors on the fit, also shown in the figure, are 
obtained by jackknifing the procedure. 

Plotted in Fig. 4(a) are the results of gg dis with the same sea-quark mass as 
those of the valence- (and cloud-) quarks in the nucleon. They suggest a non-linear 
behavior in the quark-mass. This is enhanced by our finding of a very soft form factor 
(Fig. 4(b)) which is consistent with the expectations of X-PT where the pion loop 
leads to a non-analytic behavior in m^ 2 . Furthermore, this non-linear behavior is 
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seen prominently in hadron masses when dynamical fermions are included ||22|| . For 
these reasons, we fit (N\uu + dd\N)di S with the linear plus m 1 / 2 form as for the C.I. 
and get a small x 2 ( see Fig- 4(a)). The extrapolation to the chiral limit is carried out 
in the same way as in the case of g A ||12 |. To calculate (iV|ss|./V), we fix the sea-quark 
mass at 0.154 and extrapolate the valence-quark mass to the chiral limit with the 
form C + D\Jm + m s to reflect the m 3 K dependence of the nucleon mass from the 
kaon loop in x PT. These results are also plotted in Fig. 4(a). 

From Fig. 4(a), we find that (N\uu + dd\N) dis = Z s (N\uu + dd\N)^ is = 5.41(15). 
This is 1.8(1) times the C. I. and is consistent with previous indirect calculations based 
on dM^ /dm with dynamical fermions ||, a direct calculation with staggered fermions 
||19|| , and the recent direct calculation [^] which gives a ratio of 2.2(6). Similarly, 
we find from Fig. 4(a) that (iV|ss|iV) = Z s f{ma)(N\ss\N) L = 1.53(7) where we 
have included the finite ma correction factor f(ma) = 0.79 which was computed by 
comparing the triangle diagram in the continuum and on the lattice ||18j|. This is 



much smaller than the recent calculation |20| which gives (iV|ss|./V) = 2.84(44). Part 
of the disagreement comes from the fact that a finite ma correction factor which is 
only appropriate for a C. I. was used in ref. p0[ for the D.I.. This leads to an 
overestimate by ~ 30%. In addition, summing J2x^(x,t) in eq. (Kl) over the edges 



in time where the fixed B. C. is applied as is done in Ref. |2(| gives an unphysical 
effect. Our difference might be reconciled by these two effects. 

From the above results, we list (p\uu\p), (p\dd\p), (N\ss\N), F s = ((p\uu\p) — 
(N\ss\N))/2, and D s = ((p\uu\p) - 2(p\dd\p) + (N\ss\N))/2 in Table 1. We see that 
both Ds and Fs compare favorably with the phenomenological values obtained from 
the SU(3) breaking pattern of the octect baryon masses with either linear ||, |20 



or 



quadratic mass relations p3| . Especially, we should point out that the agreement 
is significantly improved from the valence quark model which predicts Fs < 1 and 
Ds = and also those of the C. I. alone. The latter yields Fs = 0.91(13) and 
Ds = —0.28(10) which are only half of the phenomenological values || ^(J p3| . This 
underscores the importance of the sea-quark contributions. We also obtain the form 
factor gs,dis((l 2 ) = (N\uu + dd\N)% is (q 2 ) for the D.I. as plotted in Fig. 4(b). We 
see that it is exceedingly soft which is reminiscent of the two tc intermediate state 
in the x PT calculation Q. This possibility can be seen in Fig. 1(b) with two tt 
dominance. Indeed, if we assume that the D. I. part completely saturates a n N with 
g s = 8.43(24), it would give Aa^v = 11.5(2.1) MeV similar to that of the ^PT 
calculation ||]. However, there is also the C. I. part (fig. 2(b)) which is much harder 
than the D.I. When combined, it yields a scalar form factor gs{Q 2 ) which is softer 
than g\(q 2 ) and becomes close to GE^q 2 ) of the proton. They are plotted in Fig. 5 
for comparison. Fitting the gs(Q 2 ) to a dipole form gives a dipole mass mo = 0.80(4) 
GeV. This predicts Aovtv = 6.6(6) MeV, much smaller than the 15.2(4) MeV obtained 
solely based on the two-7r dominance. We conclude from this that the xPT calculation 
Q is relevant to the D.I. but missed the C.I. which maybe dominanted by the scalar 
meson. On the other hand, the (N\ss\N)(q 2 ) comes only from the D.I., hence is very 
soft. Its r.m.s. radius (?" 2 ) </ 2 (s.s) = 1.06(9) fm can be interpreted as the size of the 
KK meson cloud in the scalar channel (see Fig. 1(b)). 
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For the parameter y in eq. @, we find it to be 0.36(3). Both Aa n N and y differ 
significantly from the phenomenological solution based on xPT as mentioned earlier 
which did not take into account the C. I. with a possible scalar dominance. Our 
results on Aa n N and y strongly suggest a higher a n N = — A<j n N ~ 53 MeV, 
assuming ~ 60 MeV and a 1 - ) ~ 32 — 35 MeV. Now, we turn to our result of 
o-kN- Our direct calculation gives (N\uu + dd\N) = 8.43(24) and a n N — 

49.7(2.6) 

MeV. This is about one and half a smaller than 53 MeV inferred from Aa^N and 
y. Since the direct computation of <r n N involves the determination of the quark 
mass which is more susceptible to systematic errors (such as the extrapolation in 
the quark mass and the infinite volume limit) than the q 2 dependence of the form 
factor and the ratio y, we believe that our result on u^n is less reliable than Act^n 
and y. To examine the sensitivity of these three quantities as far as the chiral limit 
extrapolation is concerned, we fit them to a linear function in m instead of m 1//2 and 
find that Aa nN = 4.7(8) MeV, y = 0.42(3), and a wN = 39.0(2.0) MeV. Again, we 
see that both Aa nN and y favor a higher a n N ~ 55 MeV which is very close to the 
above estimate of 53 MeV with the m 1 / 2 extrapolation. Yet, the directly calculated 
ovtv falls short of this expectation and is also much smaller than that of the m 1//2 
extrapolation. 

Table 1: Scalar contents, Aa n N, J and ct^n in comparison with phenomenology 





Lattice 


Phenomenology 


(p\uu\p) 


4.55(16) 




(p\dd\p) 


3.92(16) 




(N\ss\N) 


1.53(7) 




F s 


1.51(12) 


1.52 — 1.81@ 


D s 


-0.88(28) 


-0.520 — -0.57(23| 


(r 2 )T(ud) 


0.85(4) fm 




(r 2 )T(s) 


1.06(9) fm 




Aa wN 


6.61(59) MeV 


15.2(4) MeV [|] 


y 


0.36(3) 


0.2 - 0.3 j§ 


VttN 


49.7(2.6) MeV 


45 MeV [| 


&KN 


362(13) MeV 


395 MeV 24] 



Clearly, calculations on larger lattices and smaller quark masses will be needed 
to bring the systematic errors under control and obtain a completely consistent so- 
lution on Aa n N,y, and a n N ■ Eventually dynamical fermions need to be included 
to complete the picture. Nevertheless, based on what we have on a qualitative 
and semi-quantitative level, we find that a consistent solution might be close to 
Aa n N = 6.6(6)MeV, y = 0.36(3), and a^^ ~ 53 MeV which are significantly dif- 
ferent from the present phenomenological values. We should stress that our results 
on Fs and Ds, like their counterparts in the axial couplings, agree well with those 
deduced from the SU(3) breaking pattern of the octect baryon masses and that the 
D.I. is the important ingredient for this agreement. In addition, we report the KNa 
term okn — (rh + m s )(N\uu + dd + 2ss|iV)/4 in Table 1. If we assume that they 
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are similarly depressed as in a^N, we would then predict the final ctkn at 389(14) 
MeV. It agrees with <Jkn — 395 MeV from a recent chiral analysis of KN scattering 
|24| . Finally, we note that m s (N\ss\N) = 183 (8) MeV . Together with the kinetic and 



potential energy contribution of —90 MeV p5| , the strange quark contributes about 
90 MeV to the nucleon mass. 

This work is partially supported by DOE Grant DE-FG05-84ER40154. The au- 
thors wish to thank G.E. Brown, N. Christ, T. Draper and M. Rho, for helpful 
comments. 
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Figure Captions 

Fig. 1 (a) Connected insertion, (b) Disconnected insertion. 

Fig. 2 (a) The lattice gs con for the C. I. as a function of the quark mass ma. The 
chiral limit result is indicated by •. (b) The form factor gs con {q 2 ) at the chiral limit. 

Fig. 3 The ratios of eq.(f)]) for the scalar current are plotted for the 3 Hi CclSCS. ME is 
the fitted slope. 

Fig. 4 (a) The D.I. of (N\uu + dd\N) and (A|ss|A) as a function of ma. The chiral 
limit result is indicated by •. (b) The corresponding form factors. 

Fig. 5 The normalized form factor gs(q 2 )/ds{0) compared with Ge(q 2 ) and g\(q 2 ) 
and their respective experimental results. 
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Fig. 5 
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